function [ M ] = tensorProductGradient( s, ds, A, dA )
%TENSORPRODUCTGRADIENT Summary of this function goes here
%   Detailed explanation goes here

dsA = reshape(kron(ds, A), 3, 3, 3);
dAs = dA * s;

M = dsA + dAs;

% M = dsA;

end

